home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / SPPSL.z / SPPSL
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. SSSSPPPPPPPPSSSSLLLL((((3333FFFF))))                                                            SSSSPPPPPPPPSSSSLLLL((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      SPPSL   - SPPSL solves the real symmetric positive definite system A * X
  10.      = B using the factors computed by SPPCO or SPPFA.
  11.  
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  14.       SUBROUTINE SPPSL(AP,N,B)
  15.  
  16. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  17.      On Entry
  18.  
  19.      AAAAPPPP REAL (N*(N+1)/2)
  20.         the output from SPPCO or SPPFA.
  21.  
  22.      NNNN INTEGER
  23.         the order of the matrix  A .
  24.  
  25.      BBBB REAL(N)
  26.         the right hand side vector.  On Return
  27.  
  28.      BBBB the solution vector  X . Error Condition
  29.  
  30.      AAAA division by zero will occur if the input factor contains a zero on the
  31.      diagonal.  Technically, this indicates singularity, but it is usually
  32.      caused by improper subroutine arguments.  It will not occur if the
  33.      subroutines are called correctly and  INFO .EQ. 0 .  To compute
  34.      INVERSE(A) * C  where  C  is a matrix with  P  columns
  35.         CALL SPPCO(AP,N,RCOND,Z,INFO)
  36.         IF (RCOND is too small .OR. INFO .NE. 0) GO TO ...
  37.         DO 10 J = 1, P
  38.         CALL SPPSL(AP,N,C(1,J)) 10 CONTINUE LINPACK.  This version dated
  39.      08/14/78 .  Cleve Moler, University of New Mexico, Argonne National Lab.
  40.      Subroutines and Functions BLAS SAXPY,SDOT
  41.  
  42.  
  43.  
  44.  
  45.  
  46.  
  47.  
  48.  
  49.  
  50.  
  51.  
  52.  
  53.  
  54.  
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.